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SECTION  1 


INTRODUCTION 

1.1  Background 

Investigations  of  several  recent  aircraft  accidents  indicate  that 

wind  shear  is  a factor  which  can  cause  aircraft  crashes.  One  accident 

was  a crash  during  an  ILS  approach  of  Iberia  Airlines  McDonnell  Douglas 

DC  10-30  at  Logan  International  Airport,  Boston,  on  December  17,  1973. 1 

Another  accident  was  an  Eastern  Airlines  Boeing  727  crash  on  June  24,  1975 

2 

during  an  approach  to  John  F.  Kennedy  Airport. 

During  this  time  frame,  the  "Speckled  Trout"  C-135A  flight  test 
program  at  the  Air  Force  Flight  Dynamics  Laboratory  was  investigating 
advanced  automatic  landing  techniques.  Various  inconsistencies  in 
performance  were  noted.  Wind  shear  was  suspected  as  the  cause.  Conse- 
quently, a hybrid  simulation  of  the  aircraft  and  automatic  control  system 
was  undertaken  to  investigate  the  effects  of  wind  shear  on  performance. 
This  study  indicated  that  a logarithmically  decreasing  headwind  from 
50.67  ft/sec  (30  kt)  at  510  ft  AGL  (above  ground  level)  to  0 ft/sec  at 
10  ft  AGL  resulted  in  an  automatic  landing  721  ft  short  of  the  no  wind 
case  and  a sink  rate  at  touchdown  of  -6.2  ft/sec  instead  of  -2.1  ft/sec 

3 

for  the  no  wind  case.  From  an  altitude  of  1000  ft  to  510  ft,  the  wind 
was  a constant  50.67  ft/sec.  This  particular  wind  shear  condition 
resulted  in  the  worst  performance  compared  to  steady  state  winds  and 
linear  shears  of  comparable  magnitude  (logarithmic  winds  were  considered 
more  representative  of  atmospheric  conditions).  Considering  the  above 
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logarithmic  shear,  a similar  30  kt  headwind  shear  which  is  linear,  and 
a constant  30  kt  headwind,  the  longitudinal  range  dispersion  at  touchdown 
was  2017  ft  and  the  sink  rate  dispersion  at  touchdown  was  5.55  ft/sec. 
Dispersion  is  the  span  from  the  minimum  to  the  maximum  values;  each 
extreme  limit  is  determined  by  a wind  condition.  The  longitudinal  control 
equations  were  based  on  the  conventional  operational  autopilot  installed 
in  the  aircraft;  during  the  time  based  exponential  flare,  the  throttles 
retard  linearly. 

1 . 2 Analysis  Technique 

The  objective  of  this  investigation  is  to  design  a flare  control 
law  that  will  enable  an  aircraft  on  an  automatic  approach  in  the  presence 
of  unknown  wind  shears  to  achieve  the  smallest  possible  longitudinal 
dispersion  and  sink  rate  deviation  at  touchdown  relative  to  a no  wind 
case.  The  flare  portion  of  the  automatic  landing  was  chosen  since  the 
flare  control  law  is  a major  cause  of  poor  landings  in  the  presence  of 
wind  shears."^ 

The  best  possible  performance  should  be  obtained  by  using  optimal 
control  theory.  Nonlinear  equations  of  motion  cannot  be  used  since  a 
control  law  has  to  be  obtained  for  every  new  set  of  initial  conditions 
and  the  computation  speed  is  not  great  enough  for  real  time  application. 
Linear  equations  of  motion  with  a quadratic  performance  integral  permit 
a solution  which  yields  a control  law  independent  of  initial  conditions; 
the  optimal  control  gains  are  fixed  prior  to  installation  of  an  automatic 
landing  system  in  an  aircraft  and  are  valid  for  any  flare  initial  conditions. 
In  application,  of  course,  these  initial  conditions  should  fall  within  a 
broad  set  of  appropriate  flare  initial  conditions  to  ensure  acceptable 
trajectories. 
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Linear  optimal  control  solutions  for  flare  have  been  investigated. 

However,  the  effects  of  arbitrary  wind  shears  on  longitudinal  dispersions 

and  touchdown  sink  rates  have  not  been  addressed.  Merriam4  and  Neal"*  did 

not  consider  wind  effects.  Applying  their  optimal  control  gains,  which 

are  functions  of  time,  in  a wind  environment  will  not  yield  satisfactory 

results.  Dispersions  and  sink  rates  will  change  considerably  in  the 

presence  of  winds  since  the  air  mass  is  in  motion.  As  an  example,  a 

constant  headwind  will  cause  the  aircraft  to  be  further  from  the  no  wind 

touchdown  point  at  a corresponding  point  in  time.  Huber6  included  a 

linear  wind  shear  in  his  work;  but,  the  wind  shear  was  known  a priori. 

The  wind  shear  was  placed  in  the  equations  of  motion;  optimal  control 

gains  reflected  a specific  wind.  The  solution  is  valid  for  only  that 

7 

wind  condition.  Trankle  used  an  optimal  regulator  with  integral  control 
to  track  an  altitude  profile.  This  technique,  however,  can  only  handle 
constant  winds.  The  results  are  not  satisfactory  if  wind  shears  are 
introduced. 

The  approach  in  this  investigation  was  to  use  linear  optimal  control 
theory  and  model  the  performance  index  after  a pilot's  performance. 
Specifically,  a trajectory  was  not  tracked.  Instead,  terminal  weights 
were  placed  on  pitch  rate,  airspeed,  thrust,  touchdown  point,  sink  rate, 
and  altitude.  Since  airspeed  is  critical  during  flare,  an  airspeed 
profile  was  tracked.  Also,  in  order  to  keep  normal  acceleration  small, 
in-flight  weights  were  placed  on  pitch  rate,  normal  airspeed,  and 
elevator  position.  The  controls  were  control  column  angular  displacement 
and  throttle  angular  displacement. 
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Furthermore,  a crucial  difference  between  this  investigation  and 
previous  work  that  used  time  varying  optimal  gains  is  the  application  of 
the  time  varying  gains  based  on  range  from  a no  wind  touchdown  point. 

This  scheduling  is  feasible  since  range  from  a touchdown  point  is  one  of 
the  aircraft  states.  Range  may  be  determined  in  the  future  from  a 
microwave  landing  system  (MLS)  or  global  positioning  satellite  system  (GPS). 
As  was  mentioned  previously,  gain  scheduling  based  on  time  in  the  presence 
of  wind  does  not  yield  satisfactory  results.  In  this  study,  time  varying 
gains  are  obtained  for  a no  wind  condition.  The  gains  are  stored  according 
to  the  range  of  the  aircraft  from  the  no  wind  touchdown  point.  Constant 
wind  and  wind  shears  are  then  treated  as  disturbances  to  the  system  in 
simulations  of  the  aircraft  using  the  optimal  no  wind  control  law. 

Also,  in  the  simulations,  ground  effect  is  treated  as  a disturbance. 
But,  unlike  wind  which  cannot  be  predicted,  a nominal  ground  effect  model 
is  included  in  the  solution  for  the  optimal  gains. 

A CDC  6600  digital  computer  was  used  for  all  computations.  The 
optimal  control  gains  and  simulation  are  for  a C-135A  aircraft. 

Simulation  indicated  that  the  optimal  control  gains  based  on  range 
provide  smooth  flare  trajectories  with  a longitudinal  touchdown  dispersion 
of  only  71.1  ft  and  sink  rate  dispersion  at  touchdown  of  only  .17  ft/sec. 

30  kt  constant,  linearly  shearing,  and  logarithmically  shearing  headwinds 
were  included  in  this  investigation;  these  winds  are  the  same  as  those 
mentioned  previously  in  the  discussion  of  the  "Speckled  Trout"  program. 

Thus,  compared  to  a conventional  operational  autopilot,  the  linear  optimal 
control  law  permits  a very  significant  reduction  in  longitudinal  and  sink 
rate  touchdown  dispersions  in  wind  conditions. 
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SECTION  2 
SYSTEM  MODEL 


2.1  Aircraft  Model 


Linearized  longitudinal  equations  for  perturbations  about  trimmed 
straight  and  level  flight  in  a landing  configuration  were  used  for  the 
aircraft  model.  Stability  axes  were  used:  right  hand  coordinate  system 
with  x axis  along  the  steady  state  velocity  vector  and  z axis  pointed 

g 

down.  The  steady  state  values  for  pitch  angle,  pitch  rate,  and  normal 

inertial  velocity  are  zero.  It  should  be  noted  that  even  though  steady 

state  pitch  angle  is  zero,  the  nose  is  slightly  high  in  straight  and 

level  flight  for  the  configuration  used.  The  aircraft’s  pitch  angle  is 

equal  to  perturbation  pitch  plus  the  straight  and  level  pitch  angle. 

The  actual  pitch  angle  is  needed  to  indicate  whether  or  not  the  main 

landing  gear  impacts  the  ground  prior  to  the  nose  wheel.  The  perturbation 

pitch  rate  and  perturbation  normal  inertial  velocity  are  equal  to  the 

aircraft's  pitch  rate  and  normal  inertial  velocity.  The  stability 

derivatives  for  a C-135A  aircraft  in  landing  configuration  with  a steady 

state  airspeed,  U , of  261.8  ft/sec  and  gross  weight  of  160,000  lbs  arc 

g 

listed  in  Table  1.  The  equations  of  motion  are 

6 = Mw  + M*  w + M S ♦ M u + M,„  6T  + M.  6e  + M..  H 

w as  w as  q u as  6T  fie  H 

w=Z  w +Z*w  +(U  +Z)6+Z  u + \ fie  + Z„  H 

w as  w as  o q u as  fie  H 

u=-g0  + Xw  ♦ X u + X.».  6T  + X.  fie  + Xu  H. 

w as  u as  51  fie  H 
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TABLE  1 

LONGITUDINAL  DIMENSIONAL  STABILITY  DERIVATIVES^ 


M. 

w 


6T 


<Se 


Z. 

w 


u 

Z. 

<Se 

•7 

Tf 

X 

w 

X 

u 

St 

Se 


-.6916  x 10' 
-.1466  x 10' 
-.9512 

0 

.1852  x 10" 
1.390 
-.2576 
-.708 
-.0105 
-6.84 
-.302 
9.23 
-5.113 
.0719 
-.0474 

.2013  x 10~3 
1.84 
1.812 
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It  should  be  noted  that  the  variables  on  the  left  hand  side  are  referenced 


to  an  inertial  reference  frame  and  those  on  the  right  hand  side  are 
referenced  to  a reference  frame  fixed  in  the  air  mass. 

2.2  Elevator  and  Engine  Lag  Equations 

Elevator  lag  was  modelled  as  a first  order  lag  with  time  constant 
The  equation  is 

• 1 1 

6e  = - — Se  + — 6e  . 

t x c 


The  control  parameter  is  control  column  angular  displacement,  6e  . Change 
in  elevator  angle  is  Se;  up  elevator  is  positive.  Deflection  is  relative 
to  the  stabilizer  which  moves  to  trim  the  aircraft.  Stabilizer  trim  was 


not  modelled;  its  position  remains  fixed  at  its  steady  state  value.  Steady 
state  elevator  angle  is  zero  relative  to  the  stabilizer.  Elevator  time 
constant  was  set  at  1/3  sec. 

Engine  lag  was  also  modelled  as  a first  order  lag  with  time  constant 
t The  equation  is 

ST  = - - ST  + — 6th  _ (31 , 39)  (1.5)  (750)  . 

Tth  Tth  C 


The  control  parameter  is  throttle  angular  displacement,  Sth^.  Full  throttle 
deflection  is  approximately  48  deg'.  The  last  term  on  the  right  hand  side  is 
needed  in  order  to  specify  idle  thrust  level  for  the  throttle  at 
approximately  zero  deg.  The  750  in  the  second  term  on  the  right  hand  side 
of  the  equation  is  necessary  to  get  units  of  pounds  for  the  change  in 
thrust,  6T.  ST  is  the  change  in  thrust  relative  to  the  steady  state  value 
which  is  half  of  the  total  thrust  available.  Total  thrust  is  approximately 
48,000  lbs.  The  time  constant  was  set  at  1/1.5  sec. 
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2. 3 Position  in  Space 

The  inertial  location  of  the  aircraft  is  specified  by  its  altitude 
above  the  ground  and  range  from  a specified  touchdown  point  on  a runway. 
The  center  of  gravity  of  the  aircraft  is  the  reference  point  for  these 
variables.  The  center  of  gravity  is  considered  to  be  10  feet  above  the 
ground  when  the  aircraft  is  on  the  runway.  The  equations  for  determining 
range,  r,  and  altitude,  h,  are 


and 


r = U + u 
o 


h = U 9 - w. 
o 


2.4  Winds 

The  effects  of  longitudinal  winds,  u^,  are  considered  in  this  study. 
Therefore,  the  relationships  between  inertial  and  air  mass  velocities  are 


and 


u = u + u 
as  w 


w = w 


as 

The  three  types  of  winds  considered  in  this  investigation  are  constant 
headwind,  linear  headwind  shear,  and  logarithmic  headwind  shear.  The 
constant  headwind  is  50.67  ft/sec  (30  kt)  from  1000  to  10  ft  above  the 
ground.  The  linear  headwind  shear  is  a constant  50.67  ft/sec  from  1000  to 
510  ft  above  the  ground,  at  which  point  it  starts  decreasing  linearly  to 
0 ft/sec  at  10  ft.  The  equation  for  the  linear  shear  is 
50.67  h 50.67 


Uw  = 500 


50  ' 
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The  logarithmic  headwind  shear  is  similar  to  the  linear,  except  that  the 
headwind  decreases  logarithmically.  The  equation  for  this  wind  shear  is 


50.67  In (h/510) 
Uw  ‘ " ln(51) 

These  winds  are  shown  in  Figure  1. 

2.5  Ground  Effect 


50.67. 


The  ground  effect  parameter,  H,  for  the  C-135  aircraft  is‘ 


H = .0139  exp 


120  - h 
30.6818 


This  model  is  valid  from  an  altitude  of  120  ft  to  touchdown. 

2 . 6 Full  System  Model 

Using  the  previous  equations,  the  full  system  may  be  described  by 
x = Ax  + Bu  + di  + d2 

where  the  state  vector  is 


and  the  control  vector  is 
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Altitude  Above  Ground  (ft) 


Hfr  ^ 


Figure  1.  Wind  Profiles 
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SECTION  3 


OPTIMAL  CONTROL  TECHNIQUE 

3. 1 Optimal  Control  Formulation 

The  linear  optimal  control  problem  for  the  no  wind  system  with 
terminal  and  in-flight  constraints  may  be  described  by  the  state 
equation 

x = Ax^  + Bu  + dj 

and  the  performance  index 

J = ^lx(tf)-r(tf)]TS[x(tf)-r(tf)] 

♦ hi  fUx(t)-r(t)]TQ(t)[x(t)-r(t)] 

t 

o 

♦ uT(t)R(t)u(t) }dt. 

The  desired  terminal  values,  r(t^),  are 


q(tf) 

= 0 rad/sec 

Uas(V 

= -13.8  ft/sec 

«T(tf) 

= -23,540  lbs 

r(tf) 

= 0 ft 

h(tf) 

= -2.5  ft/sec 

h(tf) 

= 10  ft. 

The  pitch  rate  constraint  minimizes  rotation  at  touchdown.  The  airspeed 
and  sink  rate  specified  are  desired  for  the  C-1.35  aircraft.  The  thrust 
constraint  specifies  idle  thrust  at  touchdown.  The  range  and  altitude 
values  identify  a desired  ground  referenced  touchdown  point. 
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The  in-flight  constraints,  r(t),  are 
q(t)  = 0 rad/sec 

w(t)  = 0 ft/sec 

u (t)  = - . 1 9 1 1 2 ft /sec 

3S 


6e(t)  = 0 rad. 

These  constraints  specify  a trajectory  with  a smooth  arresting  of  sink 
rate  and  a designated  speed  bleed  profile. 

3.2  Optimal  Control  Solution 

The  solution  to  a linear  optimal  control  problem  with  terminal  and 
in-flight  constraints  is  given  by  Huber.6  He  extends  the  optimal  linear 
tracking  technique  discussed  by  Kirk^  to  include  a disturbance  vector. 
The  control  law  is 

u(t)  = F ft)  x(t)  + v(t)  . 

Huber  discusses  in  detail  a technique  for  obtaining  F and  v. 

In  determining  F and  v,  ground  effect  was  approximated  by  using  the 
following  altitude  profile  in  order  to  completely  specify  d j : 
h = 70  exp  (- . 255t) . 

It  should  be  noted  that  d_2  was  not  included  ^n  the  computation. 

With  the  optimal  control  law  specified,  the  no  wind  state  equation 
was  integrated  by  implementing  a fourth-order  Runge-Kutta  algorithm. 

The  weights  which  resulted  in  the  best  adherence  to  the  in-flight  and 
terminal  constraints  are 
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3.3  Implementation  of  Optimal  Solution 

Based  on  the  no  wind  optimal  solution,  F(t)  and  v(t)  were  tabulated 
with  respect  to  the  range  associated  with  each  time  interval:  F(r)  and 
v(r).  Simulation  of  flare  in  wind  conditions  was  then  accomplished  with 
the  following  equations: 

x = Ax^  + Bu  + (i  i + d^ 


and 

u = F (r)  x + v (r) 

where  r is  the  state  specified  in  the  state  vector  x_.  The  disturbance 
vector  dj  includes  ground  effect  which  is  a function  of  altitude.  The 
wind  disturbance  vector,  d2,  is  also  specified  as  a function  of  altitude. 

It  should  be  noted  that  the  in-flight  and  terminal  constraints, 
along  with  their  weights,  were  not  fully  determined  prior  to  simulation 
with  winds.  Simulation  in  wind  conditions  aided  in  identifying  the 
constraints  and  weights  that  were  selected. 
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SECTION  4 


RESULTS 


4. 1 Scope 

The  optimal  flare  control  law  was  applied  in  50.67  ft/sec  constant 
headwind,  50.67  ft/sec  linear  headwind  shear,  and  50.67  ft/sec  logarithmic 
headwind  shear  environments.  Two  distinct  cases  were  simulated.  In  the 
first  case,  the  initial  conditions  for  all  the  states  were  kept  the  same 
for  all  wind  conditions;  the  initial  values  were  the  same  as  for  the  no 
wind  condition.  In  the  second  case,  the  initial  values  were  different 
for  each  wind  condition.  In  both  cases,  the  no  wind  results  were  used  as 
a reference.  The  initial  conditions  for  both  cases  are  listed  in  Table  4. 
The  initial  values  were  taken  from  a hybrid  simulation  of  a C-135  aircraft 

3 

undergoing  an  automatic  approach  to  landing  on  a 2.8  deg  glideslope.  A 
conventional  operational  autopilot  was  modelled  in  that  simulation  which 
included  glideslope  and  localizer  capture,  glideslope  and  localizer 
tracking,  and  flare.  The  values  of  the  states  at  the  beginning  of  flare 
were  used  as  the  initial  conditions.  The  same  wind  conditions  described 
in  this  investigation  were  used  to  obtain  the  initial  values  for  the 
second  case.  The  second  case,  then,  reflects  the  effect  of  winds  through 
the  whole  approach.  This  permits  a realistic  evaluation  of  the  optimal 
control  law  in  an  automatic  approach  and  landing.  The  first  case,  on  the 
other  hand,  isolates  the  performance  of  the  optimal  flare  law;  a direct 
comparison  of  the  response  to  various  wind  conditions  can  be  made.  Both 
cases  are  described  by  plots  of  the  flare  trajectories  which  are  changes 
in  altitude  versus  range,  the  nine  states  versus  time,  and  the  two  control 
variables  versus  time. 
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INITIAL  CONDITIONS  FOR  FLARE 


4.2  First  Case:  Simulation  with  Identical  Initial  Conditions 


The  flare  trajectories  for  the  first  case  are  shown  in  Figure  2. 

Flare  starts  at  70  ft.  Since  the  altitude  of  the  center  of  gravity  is 
plotted,  touchdown  occurs  at  10  ft.  The  optimal  no  wind  profile  is  a 
smooth  arresting  of  sink  rate  which  would  be  acceptable  to  a pilot.  The 
linear  and  logarithmic  headwind  profiles  are  very  similar  and  do  not  vary 
too  much  from  the  no  wind  profile.  The  steady  headwind  profile  is  several 
feet  below  the  no  wind  case.  This  is  to  be  expected  since  a path  is  not 
tracked  and  a high  wind  continues  for  the  whole  flare.  The  steady  wind 
trajectory  is  acceptable;  in  fact,  this  profile  is  similar  to  a manual 
landing  in  the  same  wind  conditions.  The  touchdown  points  for  all  winds 
are  close  to  one  another.  The  slopes  at  the  ends  of  the  trajectories 
indicate  that  the  sink  rates  are  also  similar.  More  specific  information 
about  the  profiles  and  end  conditions  is  provided  by  plots  of  the  state 
variables. 

Pitch  rate  of  the  aircraft  is  shown  in  Figure  3.  The  no  wind  curve 
stays  positive;  the  aircraft  is  rotating  up  in  the  flare  to  arrest  sink 
rate  and  decrease  airspeed.  The  optimal  computation  yields  an  oscillation 
toward  the  end  of  the  flare;  however,  pitch  rate  is  positive  and  the 
maximum  value  is  less  than  .02  rad/sec.  The  linear  headwind  curve  is  very 
similar  to  the  no  wind  curve  in  oscillations  and  overall  time  for  flare. 
This  occurs  since  the  linear  headwind  at  flare  initiation  is  not  too  great 
and  it  shears  to  zero;  its  value  is  6.08  ft/sec.  The  logarithmic  headwind 
at  flare  is  much  higher;  its  value  is  25.07  ft/sec.  This  wind  causes  a 
a considerably  slower  groundspeed  which  results  in  a longer  flare  time  in 
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order  to  meet  the  range  constraint.  The  logarithmic  headwind  also  shears 
more  quickly  in  the  flare.  Thus,  the  logarithmic  headwind  curve  differs 
noticeably  from  the  no  wind  case.  The  constant  50.67  ft/sec  headwind 
causes  an  even  slower  groundspeed  with  a correspondingly  longer  flare 
time.  Pitch  rate  toward  the  end  of  the  flare  is  smaller  than  in  the  no 
wind  case  in  order  to  keep  the  pitch  angle  of  the  aircraft  smaller;  this 
smaller  pitch  angle  permits  a stretching  of  the  path  in  time  so. that  the 
terminal  conditions  can  be  met.  In  the  above  discussion,  it  should  be 
noted  that  airspeed  is  tracked,  but  altitude  is  not  tracked. 

Perturbation  pitch  angle  is  shown  in  Figure  4.  As  indicated  by 
pitch  rate,  pitch  angle  increases  during  the  flare.  The  constant  headwind 
curve  indicates  the  stretching  of  the  flare  for  that  wind  condition;  even 
though  the  perturbation  pitch  angle  is  very  slightly  negative,  the 
positive  nonperturbation  value  yields  a positive  pitch  angle  at  touchdown 
for  the  aircraft.  A positive  value  for  pitch  angle  is  needed  to  have 
initial  ground  contact  with  the  main  landing  gear  and  not  the  nose  wheel. 

Normal  airspeed  of  the  aircraft  is  shown  in  Figure  5.  The  constant 
headwind  case  stands  out  again;  the  downward  airspeed  is  considerably 
lower  toward  the  end  of  flare  than  in  the  other  conditions.  As  before, 
this  permits  the  satisfaction  of  the  terminal  constraints. 

f I 

Perturbation  longitudinal  airspeed  is  plotted  in  Figure  6.  It  should 
be  noted  that  a perturbation  longitudinal  airspeed  profile  is  tracked. 

The  no  wind  curve  dlosely  follows  the  desired  airspeed.  The  linear  head- 
wind curve  is  again  very  near  to  the  no  wind  curve.  The  greater  shear  in 
the  logarithmic  headwind  case  is  evident;  during  the  first  few  seconds  of 
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Figure  4.  Perturbation  Pitch  Angles  with  Identical  Initial  Conditions  for  Various  Winds. 
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flare,  the  airspeed  is  about  2 ft/sec  slower  than  desired.  The  airspeed 
approaches  the  desired  profile  in  the  latter  part  of  the  flare.  However, 
lift  is  decreased  sufficiently  to  result  in  a slightly  lower  logarithmic 
headwind  trajectory  as  shown  in  Figure  2.  The  constant  wind  airspeed  curve 
is  not  very  different  from  the  desired  no  wind  curve  when  the  lower  ground- 
speed  is  considered;  it  should  be  noted  that  the  optimal  gains  are  applied 
based  on  range  from  a desired  touchdown  point.  The  terminal  airspeeds  for 
all  wind  conditions  are  near  the  desired  value  of  -13.8  ft/sec. 

Change  in  thrust  is  shown  in  Figure  7.  The  decrease  in  thrust  for 
the  no  wind  curve  is  necessary  to  attain  the  desired  airspeed  profile. 

The  increase  in  thrust  for  the  linear  and  logarithmic  headwind  curves  is 
needed  to  overcome  the  airspeed  loss  due  to  wind  shear;  the  greatest  thrust 
increase  is  for  the  larger  logarithmic  headwind  shear.  The  constant  head- 
wind curve  reflects  the  lower  groundspeed. 

Range  is  plotted  in  Figure  8.  The  spread  in  the  curves  is  due  to 
differences  in  groundspeed.  As  the  curves  indicate,  longitudinal  dispersion 
is  quite  small.  Specifically,  the  no  wind,  linear  headwind  shear, 
logarithmic  headwind  shear,  and  constant  headwind  terminal  range  values 
are  31.7,  30.1,  11.1,  and  -38.9  ft.  The  total  longitudinal  dispersion 
is  70.6  ft. 

Sink  rate  is  shown  in  Figure  9.  The  no  wind,  linear  headwind  shear, 
and  logarithmic  headwind  shear  curves  are  very  similar;  terminal  values 
are  -2.30,  -2.31,  and  -2.42  ft/sec.  The  constant  wind  curve  reflects  the 
stretching  of  the  flare  due  to  the  much  slower  groundspeed;  the  terminal 
value  is  -2.25  ft/sec.  The  total  difference  in  sink  rates  at  touchdown 
is  .17  ft/sec. 
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Figure  8.  Ranges  with  Identical  Initial  Conditions  for  Various  Winds. 


The  trajectory  profiles  with  respect  to  time  are  shown  in  Figure  10. 
The  curves  are  smooth  and  the  effect  of  decreased  groundspeed  is  evident. 
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Touchdown  occurs  at  10  ft. 

Change  in  elevator  angle  is  plotted  in  Figure  11.  The  elevator  angle 
increases  in  order  to  increase  the  pitch  angle  of  the  aircraft.  The 
elevator  angles  for  the  various  wind  conditions  are  generally  lower  than 
the  no  wind  curve;  additional  nose  up  rotation  of  the  aircraft  is  caused 
by  the  increase  in  thrust.  The  constant  wind  curve  has  a lower  maximum 
angle  in  order  to  stretch  the  flight  path.  In  order  to  meet  all  constraints, 
the  optimal  control  solution  yields  an  oscillation  in  the  elevator  toward 
the  end  of  flare.  As  the  curves  for  pitch  rate  and  perturbation  pitch 
angle  indicate,  this  oscillation  does  not  result  in  objectionable  rotation 
of  the  aircraft.  It  should  be  noted  that  the  rate  of  elevator  movement  is 
feasible. 

Aircraft  controls  are  plotted  in  Figures  12  and  13.  Change  in 
elevator  angle  lags  the  control  column  angular  displacement  shown  in 
Figure  12.  Change  in  thrust  lags  the  throttle  angular  displacement  shown 
in  Figure  13.  It  should  be  noted  that  the  jittery  behavior  of  the  controls 
which  is  evident  in  some  of  the  curves  does  not  show  up  in  the  elevator 
angle  or  thrust.  In  actual  implementation,  it  also  would  not  show  up  in 
the  cockpit  since  the  control  column  and  throttles  move  according  to  the 
output  of  servos. 

4 . 3 Second  Case:  Simulation  with  Different  Initial  Conditions 

This  case  is  a more  difficult  test  since  the  wind  shears  are  in 
effect  from  1000  ft  and  cause  a substantial  deficit  in  airspeed  at  flare 
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initiation.  The  flare  trajectories  for  the  second  case  are  shown  in 
Figure  14.  Even  with  the  airspeed  deficit,  the  optimal  control  law 
still  yields  an  acceptable  smooth  profile  that  meets  the  terminal 
constraints. 

The  pitch  rate  curves  shown  in  Figure  15  are  similar  to  those  in 
the  first  case.  The  aircraft  rotates  nose  up  in  flare  with  a slight 
pitch  rate  oscillation  near  touchdown. 

Perturbation  pitch  angle,  shown  in  Figure  16,  indicates  an  increase 
in  the  aircraft's  pitch  attitude  to  arrest  sink  rate  and  decrease  airspeed. 
The  main  landing  gear  still  impacts  prior  to  the  nose  wheel. 

Normal  airspeed  is  shown  in  Figure  17.  The  curves  are  similar  to 
the  first  case. 

Perturbation  longitudinal  airspeed  is  plotted  in  Figure  18.  The 
initial  airspeed  deficit  in  the  linear  and  logarithmic  headwind  shears 
is  substantially  decreased  after  a few  seconds  into  flare. 

The  change  in  thrust  curves,  shown  in  Figure  19,  indicate  large 
increases  in  thrust  in  order  to  overcome  the  initial  airspeed  deficit 
in  the  linear  and  logarithmic  headwind  shears.  Part  of  the  thrust 
increases,  just  as  in  the  first  case,  are  to  overcome  the  wind  shears 
during  the  flare. 

Range  is  shown  in  Figure  20.  The  no  wind,  linear  headwind  shear, 
logarithmic  headwind  shear,  and  constant  headwind  terminal  range  values 
are  31.7,  30.0,  12.5,  and  -39.4  ft.  The  total  longitudinal  dispersion 
is  71.1  ft. 
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Figure  16.  Perturbation  Pitch  Angles  with  Different  Initial  Conditions 
for  Various  Winds. 
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Figure  20.  Ranges  with  Different  Initial  Conditions  for  Various  Winds 


Sink  rate  is  plotted  in  Figure  21.  The  no  wind,  linear  headwind 
shear,  logarithmic  headwind  shear,  and  constant  headwind  terminal  sink 
rates  are  -2.30,  -2.26,  -2.42,  and  -2.25  ft/sec.  The  total  difference 
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in  sink  rates  at  touchdown  is  .17  ft/sec. 

The  trajectory  profiles  with  respect  to  time  are  shown  in  Figure  22. 
The  curves  are  smooth  with  touchdown  occurring  at  10  ft. 

Change  in  elevator  angle  is  shown  in  Figure  23.  The  curves  are 
similar  to  the  first  case. 

Aircraft  controls  are  plotted  in  Figures  24  and  25.  In  the  linear 
and  logarithmic  headwind  shears  in  Figure  25,  throttles  increase  to 
maximum  in  the  first  part  of  flare  in  order  to  overcome  the  initial 
airspeed  deficit. 


Figure  21.  Sink  Rates  with  Different  Initial  Conditions  for  Various  Winds 


Figure  23.  Changes  in  Elevator  Angle  with  Different  Initial  Conditions 
for  Various  Winds. 


0-2i  - 


24.  Control  Column  Angular  Displacements  with  Different  Initial 
Conditions  for  Various  Winds. 


Figure  25.  Throttle  Angular  Displacements  with  Different  Initial  Conditions 
for  Various  Winds. 


SECTION  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  Conclusions 

Initially,  considering  the  no  wind  environment,  the  linear  optimal 
control  technique  presented  in  this  investigation  yields  a smooth  flare 
to  touchdown.  To  achieve  this  flare,  in-flight  and  terminal  constraints 
were  placed  on  the  aircraft.  The  in-flight  constraints  are  on  pitch  rate, 
normal  airspeed,  perturbation  longitudinal  airspeed,  and  change  in  elevator 
angle.  The  terminal  constraints  are  on  pitch  rate,  perturbation 
longitudinal  airspeed,  change  in  thrust,  range,  sink  rate,  and  altitude. 

Furthermore,  and  most  importantly,  the  optimal  control  law  developed 
in  the  no  wind  environment  will  yield  smooth  flares  to  touchdown  with  very 
small  terminal  range  and  sink  rate  dispersions  in  the  presence  of  large 
constant  headwinds,  linear  headwind  shears,  and  logarithmic  headwind  shears. 
The  optimal  time  based  gains  obtained  in  the  no  wind  environment  are  applied 
with  respect  to  range;  there  is  a one  to  one  correspondence  between  range 
and  time.  The  headwinds  applied  as  disturbances  were  50.67  ft/sec  from 
1000  to  510  ft  above  the  ground;  from  510  to  10  ft  the  winds  stayed 
constant,  decreased  linearly  to  0 ft/sec,  or  decreased  logarithmically  to 
0 ft/sec  depending  on  the  wind  being  simulated.  The  touchdown  range 
dispersion  is  only  71.1  ft;  the  sink  rate  dispersion  is  only  .17  ft/sec. 

As  a comparison,  a simulation  of  a conventional  operational  autopilot  in 
the  same  wind  conditions  resulted  in  a longitudinal  dispersion  of  2017  ft 
and  sink  rate  dispersion  of  5.55  ft/sec  with  a maximum  sink  rate  of 
-6.2  ft/sec. 
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The  optimal  linear  control  gains  can  be  readily  implemented  in  actual 
flight  operations  since  only  a few  multiplication  calculations  are  required 
to  obtain  the  control  signal  and  the  gains  are  good  for  any  initial 
conditions.  It  should  be  noted  that  nonlinear  optimal  control  cannot  be 
implemented  due  to  the  long  computation  time  required  to  obtain  a solution 
and  the  need  to  recompute  gains  for  each  set  of  initial  conditions. 

5.2  Recommendations 

Since  airspeed  change  during  flare  is  substantial,  the  effect  of 
changes  in  stability  derivatives  should  be  investigated.  Also,  a nonlinear 
simulation  of  the  aircraft  with  linear  optimal  control  gains  should  be 
accomplished  prior  to  flight  testing  in  order  to  verify  all  of  the 
linearizing  assumptions.  This  simulation  should  include  a method  to  handle 
stabilizer  trim.  Also,  tailwinds  should  be  introduced  to  check  the 
aircraft's  response  in  these  winds;  the  throttles  should  come  off  in 
tailwind  shears  since  airspeed  increases. 

The  response  of  the  aircraft  to  gusts  should  be  studied.  It  should 
be  noted  that  a function  of  normal  acceleration  is  weighted  in-flight. 

Optimal  gains  for  various  aircraft  configurations,  weights,  and 
airspeeds  will  have  to  be  obtained  prior  to  flight  test.  In  this  study, 
only  one  aircraft  condition  was  used.  Provisions  for  efficiently  storing 
the  various  sets  of  gains  in  a digital  flight  computer  should  be 
investigated. 

As  this  study  indicates,  change  of  thrust  is  extremely  important  for 
small  range  and  sink  rate  dispersions  when  in  the  presence  of  wind  shears. 
In  order  to  get  a satisfactory  comparison  to  classical  flare  techniques, 
the  classical  method  should  employ  a very  good  automatic  throttle  system. 
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In  order  to  have  a smaller  airspeed  deficit  at  flare  initiation 
when  winds  are  present,  improvements  in  glideslope  tracking  should  be 
investigated.  Also,  localizer  tracking  and  alignment  should  be 
investigated  from  the  standpoint  of  effects  to  the  optimal  flare. 
Specifically,  the  nonlinear  simulation  that  was  mentioned  previously 
should  also  include  the  lateral  equations  with  crosswind  disturbances. 

It  should  be  noted  that  the  problem  that  was  solved  in  this  study 
is  not  just  an  improvement  in  aircraft  flare  performance.  The  optimal 
control  technique  presented  here  is  good  for  placing  an  air  vehicle  at 


any  point  in  space,  with  any  terminal  conditions,  with  any  in-flight 
constraints,  in  any  deterministic  wind  conditions;  this  can  be  done 
regardless  of  the  vehicle's  initial  conditions. 
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